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LONG-TERM  GOALS 

To  provide  a  capability  to  extend  limited  observations  of  local  waves  and  currents  into  a  larger  domain 
using  physics-based  models  within  the  littoral  zone.  It  is  anticipated  that  UUVs  and/or  UAVs  would 
be  used  to  provide  the  local  observations,  and  the  tool  would  be  able  to  be  run  quickly  on  a  laptop. 

OBJECTIVES 

The  objectives  of  the  proposed  program  are  to  utilize  video  from  UAVs  in  the  littoral  region  to 
generate  estimates  of  wave  conditions  and  bathymetry.  This  would  be  the  first  step  to  the  more  long¬ 
term  goals,  which  would  require  incorporating  UUV  and  in  situ  observations  as  well  as  estimating 
currents.  The  proposed  program  will  have  three  parts. 

(1)  Develop  a  ray-trace  based  solver  for  the  wave  action  balance  equation  that  will  allow  prediction  of 
shoaling  wave  spectra  given  the  local  bathymetry  and  currents,  and  will  allow  the  assimilation  of  wave 
observations  inside  the  domain  into  the  wave  field  coming  into  the  domain  that  would  reproduce  the 
observations. 

(2)  Using  the  ray-trace  solver  developed  above,  generate  estimates  of  wave  conditions  and  bathymetry 
from  video  observations  of  breaking  waves. 

(3)  Using  the  ray-trace  solver  developed  above,  generate  estimates  of  wave  conditions  from  video 
observations  of  non-breaking  waves  (i.e.  seaward  of  the  surf  zone). 

APPROACH 

(1)  Development  of  the  ray-trace  wave  action  equation  solver 

Previous  ONR  programs  have  developed  solvers  for  the  wave  action  balance  equation,  which  describes 
the  evolution  of  wave  spectra  as  the  wave  propagate  through  a  current  field  and  through  local 
bathymetry  changes.  For  this  program,  these  codes  will  be  modified  to  include  shoal  conditions,  and 
the  energy  balances  in  the  code  will  be  able  to  re-produce  what  the  SWAN  model  has  as  sources  and 
sinks  of  wave  energy.  The  assimilation  code  will  be  based  on  standard  variational  calculus  to  develop 
the  gradient  of  the  error  (mis-match  between  simulated  and  observed  wave  conditions)  with  respect  to 
the  wave  field  at  the  boundary,  and  then  a  conjugate  gradient  search  engine  will  be  developed  that 
utilizes  this  gradient.  This  will  be  tested  using  simulated  data  and,  if  available,  in  situ  wave 
observations  of  shoaling  waves  from  previous  ONR  experiments  (including  Sandy  Duck  '97  and 
NCEX). 
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(2)  Assimilation  of  video  observations  of  breaking  waves 

The  wave  aetion  equation  solver  will  eontain  loss  of  energy  due  to  breaking  during  shoaling.  A 
forward  model  for  optieal  intensity  as  a  funetion  of  breaking  will  be  developed,  and  then  used  to  find 
the  loeal  bathymetry  and  wave  height  that  ean  reproduee  the  observed  breaking.  We  antieipate 
working  elosely  with  NRL/SSC  to  utilize  their  pole-based  video  data  of  breaking  waves  in  a  number  of 
loeations  as  test  data  for  validation.  In  addition  NRL/SSC  will  provide  UAV  video  data  as  an 
additional  test  set. 

(3)  Assimilation  of  video  observations  of  non-breaking  shoaling  waves 

A  forward  model  for  generating  optieal  modulations  for  non-breaking  waves  will  be  ineluded  into  the 
system.  This  will  be  used  to  assimilate  video  of  waves  seaward  of  the  surf  zone  into  estimates  of  wave 
heights.  NRL/SSC  is  developing  a  bathymetry  algorithm  based  on  wave  dispersion  that  we  antieipate 
ineorporating  and  merging  with  the  surf  zone  algorithms  developed  earlier.  The  final  eode  will  be  able 
to  ineorporate  UAV  observations  throughout  the  littoral  region. 

WORK  COMPLETED 

This  program  started  late  in  CY  '07,  so  the  first  year  is  not  yet  eomplete.  We  have  developed  the  ray- 
traee  solver  to  the  wave  aetion  balanee  equation  for  shoaling  waves  and  ineorporated  the  sourees  and 
sink  terms  that  are  in  the  SWAN  model.  This  required  some  modifieations  to  the  SWAN  approaeh, 
sinee  the  ray-traee  solver  does  not  have  the  entire  wave  speetrum  available  at  eaeh  time.  This  means 
that  we  need  to  have  a  different  method  of  ineorporating  the  speetral  integral  parameters  that  are  in 
SWAN.  We  use  sealed  versions  of  the  equilibrium  speetrum  integrals  (whieh  we  ealeulated  initially) 
where  the  seale  faetor  eomes  from  a  eomparison  of  the  speetrum  value  for  the  given  speetral  loeation 
we  are  ealeulating  versus  the  equilibrium  speetrum  value  at  the  same  frequeney  loeation.  We  have 
done  a  eomparison  between  SWAN  and  this  model  for  a  linear  beaeh  and  a  quasi-monoehromatie 
wave  and  have  shown  that  they  are  identieal  (i.e.  we  ean  reprodueed  the  SWAN  result).  We  are 
eurrently  performing  a  eomparison  between  SWAN  and  our  model  for  a  more  realistie  beaeh  (the 
aetual  bathymetry  measured  at  Duek,  N.C.)  and  a  more  realistie  wave  speetrum  (a  Pierson-Moskowitz 
wind-wave  speetrum).  We  are  also  in  the  proeess  of  eontaeting  other  ONR  investigators  to  see  if 
historieal  data  of  wave  speetra  from  shoaling  waves  are  available  from  previous  experiments  as 
validation  oases.  In  partioular  we  are  foousing  on  Sandy  Duek  '97  and  NCEX.  If  available,  we  will 
use  these  data  to  oompare  shoaling  wave  speetra  to  model  outputs.  Before  the  end  of  the  first  year  we 
will  also  develop  the  assimilation  eode  that  will  take  observations  of  wave  speetra  within  the  domain 
and  determine  the  inooming  wave  field  that 

RESULTS 

We  have  oompleted  the  ray-traee  solver  for  the  wave  aetion  balanee  equation  and  oompared  results  to 
SWAN  for  a  linear  beaeh  and  a  quasi-monoehromatie  wave  field.  We  generated  wave  speetra  as  a 
funetion  of  distanee  from  shore  (and  thus  water  depth)  for  eaeh  model  and  ealeulated  the  peak  wave 
length  (via  a  weighted  integration  of  the  speetrum)  from  eaeh  set  of  speetra.  A  eomparison  is  shown  in 
Figure  1  where  there  is  essentially  no  differenee  between  SWAN  (blue  line)  and  the  ray-traee  eode 
(referred  to  as  Wabe  in  the  figure).  Figure  2  then  shows  a  eomparison  of  signifieant  wave  height 
between  the  two  models.  The  ray-traee  eode  (again  referred  to  as  Wabe  in  the  figure)  was  run  with 


2 


two  different  values  for  the  gamma  parameter  that  is  in  the  SWAN  energy  dissipation  model  for 
breaking  waves.  Note  that  the  value  of  0.6  essentially  reproduces  the  SWAN  plot.  Finally,  Figures  3 
and  4  shows  a  comparison  of  the  spectra  locations  and  shapes  for  different  water  depths.  Although  the 
scaling  of  the  spectra  are  different  (due  to  the  way  they  are  output)  the  locations  of  peaks  in  wave 
number  and  the  shape  of  the  spectra  are  essentially  identical. 


Figure  1:  Comparison  of  the  peak  wave  length  of  the  spectra  as  a  function  of  the  water  depth  for  the 
spectra.  The  SWAN  model  results  are  in  blue  and  the  ray-trace  model  results  (referred  to  as  Wabe) 
are  in  pink  and  green.  Note  that  they  are  essentially  equivalent. 


Figure  2:  Comparison  of  significant  wave  height  (Hmo)  between  SWAN  (blue  line)  and  the  ray- 
trace  code  (Wabe)  for  two  different  values  of  the  gamma  parameter  in  the  SWAN  energy  dissipation 
model  for  breaking  waves.  Note  that  the  Wabe  run  with  gamma  =  0.6  is 
essentially  identical  to  SWAN 
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Figure  3:  SWAN  model  output  showing  cuts  through  the  wave  spectra  for  various  water  depths  (and 
thus  various  distances  from  shore).  Depths  of  5m  and  deeper  are  essentially  the  same,  then  the 
spectrum  starts  to  shift  in  wavenumber  and  lose  energy  for  2.6m  and  1.1m  depths. 


Figure  4:  Same  as  Figure  2  but  for  the  Wabe  code.  Note  however  that  the  spectra  are  not  scaled 
the  same,  this  is  strictly  to  validate  the  shape  of  the  spectra. 
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As  mentioned  above,  we  are  eurrently  performing  the  same  comparison  for  more  realistic  waves  and 
bathymetry. 

IMPACT/APPLICATIONS 

If  successful,  these  tools  will  be  used  by  the  warfighter  to  provide  a  wide-area  view  of  current  littoral 
wave  conditions  during  operations,  thus  increasing  the  field-of-view  of  local  assets  (such  as  UAVs). 

RELATED  PROJECTS 

There  are  no  ongoing  related  projects  that  are  closely  identified  with  this  project. 
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